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Summary 


An analytical study has been conducted to predict the effect of an oscillating 
stream on the time dependent sectional pressure and lift coefficients of a model 
propeller blade. The assumption is that as the blade sections encounter a wake, the 
actual angles of attack vary in a sinusoidal manner through the wake, thus each 
blade is exposed to an unsteady stream oscillating about a mean value at a certain 
reduced frequency. On the other hand, an isolated propeller at some angle of attack 
can experience periodic changes in the value of the flow angle causing unsteady 
loads on the blades. Such a flow condition requires the inclusion of new expressions 
in the formulation of the unsteady potential flow around the blade sections. These 
expressions account for time variation of angle of attack and total shed vortices in 
the wake of each airfoil section. It was found that the final expressions for the 
unsteady pressure distribution on each blade section are periodic and that the 
unsteady circulation and lift coefficients exhibit a hysteresis loop. 


Introduction 

Propeller performance 1 and radiated noise^ can be substantially altered by 
the exposure to a non-uniform freestream, such as the flow produced by the wake of 
an upstream wing or pylon. It has also been shown that an isolated propeller at some 
angle of attack is also exposed to the unsteady loads experienced by the propeller 
blade sections^. While many researchers, motivated by unsteady flow encountered 
on helicopter rotor blades or inside turbomachinery, have studied the physics of 
airfoils exposed to nonuniform flow^-S, the present analysis addresses the problem of 



propeller aerodynamics from a technology base point of view. This analysis is 
centered on the application of ANOPP6 (Aircraft Noise Prediction Program), 
developed at NASA Langley Research Center, to the unsteady flow environment. 

While ANOPP is capable of coupling the aerodynamic calculations to the analytical 
noise prediction segment, only the aerodynamic computations leading to the 
unsteady blade section pressure distribution and their unsteady lift coefficients are 
considered here. The derevation of the equations for unsteady terms are descibed in 
the analysis section of the report and in the final portion of the analysis, closed form 
solutions have been presented for the special case of a thin airfoil. The thin airfoil 
approximation has been used in a hybrid propeller performance code developed by 
the present author. The above hybrid code has shown to be capable of computing the 
performance characteristics of advanced propellers in a very efficient manner with 
good accuracy 7. It is intended to extend the capabilities of the hybrid code to the 
unsteady environment by the use of the present two dimensional analysis in a future 
work. 


Symbols 

A terms due to complex velocity potential in Z-plane 

a radius of the mapping circle, m/s (ft/s) 

B terms due to complex velocity potential in the £-plane 

b Joukowski transformation parameter, re a 

c chord length, re a 
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section lift coefficient 

section pressure coefficient 

complex velocity potential in z-plane, re aVoo 


Mach number 


transformation. 


dC_ C 

dC’ dz 


nondimensional time, 

Uoo 

free stream velocity as seen by blade section, m/s (ft/s) 


complex velocity, m/s 
vortex location in Z-plane, re a 
complex plane of the airfoil, re a 
section angle of attack, deg. 
unsteady circulation, re aV,*, 

strength of the vortex per unit length of the wake 
increment of angle of attack 
0 - t,, rad. 

angular measure in near circle, plane 
perfect circle, £ = b eV 0 + i6 


3 


near circle, = b eV+27CiJ; 


c 

0 angular measure in £-plane 

p location of vortex in the £-plane 

<|> nondimensional potential, re aVoo 

CO reduced frequency, re — 

a 

£2 propeller angular velocity, rad./s 

Subscripts: 

le leading edge 

te trailing edge 

0 steady term 

1 unsteady components multiplying Aa coscot 

2 unsteady components multiplying Aa sincot 


Analysis 

The derivation of expressions for time dependent section pressure and 
coefficients are described in the following using the potential flow analysis. 


lift 

First, 
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expressions are described for a generic airfoil; next, special cases, such as a family of 
Theodorson's airfoils and in particular flat plate approximations, are discussed. 


a) Unsteady Potential Flow 

It has been shownl, that when the blade sections of a propeller encounter a 
wake, the blade section angle of attack can be approximated to vary in a sinusoidal 
manner throughout the passage. Thus, each blade is exposed to an unsteady 
oncoming flow. This flow can be taken as a uniform stream oscillating about a mean 
value at a certain reduced frequency, 

a = ocq + Aa sincot (1) 

Such a flow condition requires inclusion of new expressions in the 
formulation of the unsteady potential flow around the blade sections. These 
expressions are needed to properly account for time variation of angle of attack and 
total shed vortices in the wake of each airfoil section. Therefore, the effect of the 
vortical wake behind each section is included in terms of the bound circulation 
which gives rise to the magnitude of the total lift. 

It is, therefore, possible to extend the steady analysis to the problem of an 
oscillating flow about an airfoil. The unsteady potential flow is formulated for the 
perfect circle in the £ -plane. Theodorsen’s transformation is utilized to obtain the 

unsteady complex velocity in the Z-plane (see Fig. 1). The wake is modeled by 
distributed vortices composed of pairs of an isolated vortex and its image vortex 
outside and inside of the unit circle, respectively. Thus the complex velocity 
potential in the ^-plane is written: 
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Ffc) = V 00 |ce' ia +^-e i a J 




a 



Xt 6 


( 2 ) 


where the above dimensional quantities can be normalized with respect to the mean 
free stream velocity, V^, and radius of the circle, a: 


r = r/aV^ , x = x/a, y = y/Va* and £ = ae‘ 0 


Therefore, the complex velocity in the £-plane becomes 


, n w(C) 

W (g=-TT— = ie- 


i 9 


2sin(9 


a)+ £ + 


1 Q* 


Y dx 


(3) 


Xte 


with 


Q = 


2(1 - pcos0) 

2 

(1+p - 2pcos0) 


(4) 


Assuming that the distributed vortices all lie on the x-axis extending to infinity, the 
unsteady pressure coefficient can be written as 



( 5 ) 
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with 


|w(z )| 2 = w(z) * w(z), 


( 6 ) 


dC dz 


( 7 ) 


and 


w(z) = w(Q • T 

Note that <() is the real part of the complex potential, F(£). In the above expression C,' 
is the plane of the near circle obtained by the Joukowski transformation: 


z=C' + 


c 


( 9 ) 


Since the instantaneous value of angle-of-attack varies in a sinusoidal manner, it is 
logical to assume that the circulation can be expressed in the following form: 

r[t) = r 0 + AaFj coscot + Aar 2 sincot (10) 


This expression allows for any phase lead or lag in the final expression for 
circulation and, hence, lift value. Therefore, this is an adequate expression to 
describe the vorticity since it can be understood as the first term in a series 
representation for T . Since for the range of velocity defects encountered, the values 
of A a remain very small, terms of order of (Acx)^ or higher can be neglected. With 
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this in mind, it can be shown that for trailing vortices to be left behind, as the airfoil 
moves forward, the strength of the distributed vortices can be related to the bound 
circulation so that: 


y = Yi Aa cos cot + 72 Aa sin cot 


( 11 ) 


and 




x - x 



( 12 ) 

(13) 


Obviously, an expression can be found for w(£) representing a steady term and the 
unsteady contributions by using Eq. (11) in Eq. (3) to obtain 


w(Q 


= B 0 + B! Aa coscot + B 2 Aa sincot 


(14) 


where: 


B n = ie 


. -ie 


p 
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(16) 


8 



(17) 


Bo= ie 


ie 


v 

— - 2cos(e - aj + to [r l • cu - r 2 * Su) 
2k 


and terms Su and Cu are given as: 


Su= J sincojx - x t e )Q(p,0) dx 



(18) 


Cu = J cos© (x - x t ^Q(p,0) dx (19) 


Similarly, the following integrals are defined as: 


CNu = J coso) (x - x t e ) v (p,e) dx (20) 

Xtc 


SNu = | sino)(x - x t e ) v(p,0)dx (21) 



where V is the argument of the expression: 


W L = | 


( 22 ) 
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= e1e%)/(e*-I) 


These are used in the expression for time rate of change of potential: 


— = <|>' x Aa costot + <J>' 2 Aa sincot 
Bt 


C p = <t>’i = to {2(0-060) + ^ (e-ao) - colrxCNu-^SNu)) 


c|)' 2 = co / - co ( r j SNu - r /: Nu ) 


With eq. (15), (16), (17) and (7) let: 

A 0 = B 0 T 
Al=B!T 

a 2 = b 2 t 

Then the final expression for the unsteady pressure distribution is obtained as 
steady term plus unsteady components in the usual manner 

C P = C P0 + AaCpj coscot + AaCp 2 sincot 


( 23 ) 

( 24 ) 

( 25 ) 

(26) 

(27) 

(28) 
a 

( 29 ) 
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with 


C po = 1 - A 0 Ao 

(30) 

C pi = -2<t>' j — Aq Aj- AoAj 

(31) 

C p2 = -2<j>' 2 - A 0 A 2 - Ao A 2 

(32) 


where A stands for the conjugate of A. 

In the above expressions, values of the unsteady circulation components Tq, 

r j , T2 need to be specified before the pressure coefficients can be computed. The 
correct amount of the circulation is, however, related to the vorticity production and 
dissipation which itself is inherently a viscous phenomenon. For steady flow past 
thin airfoils at small angles of attack, the Kutta condition is merely a specification of 
the points of separation at the trailing edge. For unsteady flows, the rate of change 
of bound circulation must be related to the vorticity flux into the wake at the points 
of flow separation (on both upper and lower surfaces). This relation is known as the 
MRS (Moore, Rott, Sears) Criterion. This procedure would require computation on the 
viscous boundary layers on both surfaces up to the points of separation^. As was 
indicated above for airfoils at moderate angle of attack, the Kutta condition can be 
understood as trailing edge separation where the removal of a singularity due to 
transformation specifies the circulation value uniquely. Therefore, an unsteady 
extension of the Kutta condition can be understood as trailing edge separation for all 
times. This requirement specifies the circulation components uniquely, when w(£) is 

forced to vanish at a point corresponding to the trailing edge. 

r 0 = 47t sin (ct 0 - (33) 
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( 34 ) 


P 2o) Cut cos (Ote - «oj 
1 (1/27C - to Sut) 2 +(o) Cut} 2 

2(1/271 - co Sut) cos^Sjg— (Xq) 
(1/271 - co Sut) 2 + (co Cut) 2 


(35) 


In this manner then, the components of the unsteady circulation can be used to 

obtain a final expression for the unsteady lift coefficient as an unsteady extension of 

Blasius relation: 

C/ = Q 0 + AaC/j cos cot + AaC/ 2 sin cot (36) 

Indeed, the above expression will account for any phase lead or lag in the final 
expression for total lift on each blade section, to order of Aa, as the sections 

experience the oscillating streams which are caused by blade wake encounter, 

b) Family of Theordorsen Airfoils 

In the previous section, the expressions for lift and pressure coefficients were 
obtained in terms of the mapped flow past the unit circle and the metrics of the 
mapping function. In numerical analysis these constitute the Jacobian of the 
transformation. In analytical methods, such as for a family of Theodorsen's airfoils, 
the airfoil is mapped into a near circle and, subsequently, the near circle is 

mapped into the unit circle. Thus, for a known airfoil section, the coordinates and 
leading edge radius of the airfoil are given and the task is to obtain parameters of the 
near circle and the radius of the circle in terms of the airfoil data. 
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It can be shown (ref. 9) that if 


z 



for the near circle: 


C = b 

and the perfect circle is described by: 

£ = beV°e i0 = ae i0 


(37) 


(38) 

(39) 


One can write the deviation of the near circle from the perfect circle as. 

e = (K ( 4 °) 

In order to compute the deviation, £, and radius of the perfect circle, a, in terms of z 
x + iy, the relation between £ and is obtained by division of equation (38) by (39): 

£' = £ e (¥“¥o) e i(^-0) (4 1 ) 


This can be expanded in a Fourier Expansion as: 
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C’ =ex P 


X (A„ + B n ) / a n e‘ n 6 

Vn= 1 


(42) 


Thus, coefficient of the expansion can be written as: 
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^. = 1 1 \j/ cos n0 d0 


— jf- = — | \j/ sin n0 d0 


(43) 
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In an iterative manner, then, real and imaginary parts of (42) are computed 
for the following relations: 


V- Vo 


^ cos n0 + 5 jl sin n0 ) 
a n a n 


(46) 


oo 



In the limiting case for 0 = K, the leading edge radius is obtained as: 


r le = 


2 b sinh 2 \|r 
cosh \jr 


(48) 


and the Joukowski transformation parameter per chord length becomes: 

E-l.Ss. (49) 

4 8 


and the expression for the position of the isolated vortex, per radius of unit circle, 
becomes: 


(x + fx 2 - 4b 2 ) 

2 


Note that in the above expression: 


b = b * c 


(50) 


(51) 
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c) Flat Plate Approximation 


In many cases, the airfoil sections for advanced turboprops are very thin 
airfoils and the calculations above can be approximated for a flat plate airfoil section. 
In this case, the limiting values for influence functions becomes: 


Q= 


_L_ 

1-P 


(52) 


o r 


Let r\ = x-2 then: 


Q= 


-4 

(x-2 + V x 2 - 4) 


(53) 


Cu 0 = -4 1 


cos (wn) 


r\ + Vri 2 + 4r| 


dr) 


(54) 


Su 0 = -4 J 

o 


sin (ojt|) 


Tl + 


Vt | 2 + 4t) 


dri 


(55) 


Thus, the expression for the unsteady lift coefficient becomes: 

Q (t) = T(t)/2 (56) 

and expressions for C/ , Q , and are written as: 

Q = 27t sin (Xq (57) 
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(58) 


Q 


CO Clip COS OCq 

(1/2 7i - to Suq) 2 + (g)Cu 0 J 2 


Q 


2 


(1/2 7i - (0 S u 0 ) cos a D 
(1/2 n - co SuoJ 2 + (toCuo) 2 


(59) 


Results and Discussions 

The analysis described above has been employed to predict the steady state and 
unsteady lift coefficient components for blade sections of a straight bladed scale 
model (SR-2) propeller. The physical dimensions and the operating conditions are 
given in reference 1, and they are not repeated here. However, the practical 
application of the above analysis is shown by the variation of the unsteady lift 
coefficient for a value of mean angle-of-attack as described by reference 1. 
Accordingly, the section angle of attack varies in a sinusoidal form due to the passage 
of the propeller blade through the wake of the pylon. The pylon, in this case, is 
positioned upstream of the propeller in a so-called pusher configuration. The wake 
momentum deficit can cause an immediate change in angle of attack as it is seen by 
the propeller blade section. This change is periodic and can be approximated by 
sinusoidal behavior. Thus, the amplitude and reduced frequency of the motion are 
obtained from the kinematics of blade/wake interaction. A typical variation of angle 
of attack is illustrated in Fig. 2. It can be seen that the amplitude of oscillation is of 
order of 1*. Thus, the present analysis, which is restricted to small Aa, is valid for 

this type of application. 
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Next, consider the unsteady lift behavior shown in Fig. 3 for a blade section at 
37% radial distance from the propeller pitch axis. The lift coefficient was obtained 
by conformal mapping of Theodorsen’s method. A modification of the aerodynamic 
module of ANOPP was required to obtain the unsteady components of lift, drag and 
pressure coefficients under an unsteady environment. Note that the hysteresis 
behavior of the unsteady lift coefficient is recovered by this analysis without 
inclusion of unsteady boundary layer analysis. This can be explained, as discussed 
earlier, by the fact that for airfoils at small angles of attack, the unsteady Kutta 
condition is understood as trailing edge separation and uniquely defines the 
magnitude of circulation and consequently the lift curve behavior. This has 
paramount implications in understanding structural and acoustic problems. 

Reference 2 points out that lift hysteresis of aerodynamic loads, obtained under 
nonuniform inflow conditions, increases loading noise and the associated radiated 
noise by a propeller blade. 

Consider the results for the special case of a flat plate airfoil, described above 
under part c) of the analysis. This case is directly applicable to the work of reference 
7. In this work, the present author has developed a code to predict the performance 
of advanced propellers by a simple strip method, where the induced blade velocities 
are obtained by an iterative procedure. Thus, sectional angle of attack and Mach 
numbers are obtained through direct computation of section aerodynamic 
characteristics using thin airfoil theory. This analysis has shown excellent 
agreement with wind tunnel results for advanced propellers, such as SR-7. In this 
case, the contribution of the thin airfoil to the unsteady lift coefficient can be 
approximated by equations 58-60. A parametric behavior of the unsteady lift 

coefficient is simulated in Fig. 4. The angle of attack was chosen to be the same as 
shown in Fig. 2. The values of unsteady lift coefficient, Q are plotted for various 
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reduced frequencies, CO. It can be seen that the lift hysteresis loop becomes smaller 

with increasing value of w and changes the inclination of the loop accordingly. In 

addition. Fig. 5 shows that for higher mean angle-of-attack the lift coefficient and 

the size of the hysteresis loop both decrease with (Xq relative to their mean value. 

Furthermore, there is no apparent change in the direction of the loop. This is in 
accord with the explanations of reference 8, since the present analysis is limited by 
the potential flow assumption and cannot account for flow separation. consequently, 
the above results do not apply to airfoils near stall. Nevertheless, the analysis 
captures the essence of excess loads encountered under time dependent flow 
conditions. 


Conclusions 

A method of solutionhas been developed to compute the unsteady potential flow 
past airfoil sections of arbitrary shape. Even though the study is limited to 
perturbations of angle of attack from a mean value, the analysis accounts for any 
lead or lag, of the order of small amplitudes, in the final expressions for the unsteady 
aerodynamic coefficients. These coefficients then properly represent the hysteresis 
behavior encountered in a periodic flow such as the flow over advanced propeller 
blades in an unsteady flow environment. 
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